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Abstract 

Extending investigations of Antman&Malek-Madani, Schecterfe Shearer, Slemrod, 
Barker&Lewicka&Zumbrun, and others, we investigate viscoelasticity models with sur- 
face energy effects. We prove the existence of nonconstant planar periodic standing 
waves in viscoelasticity models with surface energy terms by variational methods and 
phase-plane analysis, for deformations of arbitrary dimension. For our variational analy- 
sis, we require that the mean vector of the unknowns over one period be in the elliptic re- 
gion with respect to the corresponding inviscid (i.e., elastic) model. For our (1-D) phase- 
plane analysis, we have no such restriction, obtaining essentially complete information 
on the existence of nonconstant periodic waves and bounding homoclinic/heteroclinic 
waves. Our variational framework has implications also for time-evolutionary stabil- 
ity, through the link between the action potential for the traveling-wave ODE and the 
relative entropy for the time-evolutionary system. 

Keywords: viscoelasticity, surface energy, periodic wave, Hamiltonian system. 

1 Introduction 

Viscoelasticity models are models taking both the viscosity and elasticity effects of the 
materials into account when we describe the physical phenomena in mathematical modeling. 
Since they have a wide range of applications, the study of viscoelasticity has been an 
important topic (see [AB, AM, BLeZ, D, FP, K, RZ, SI, S2, S3, SS, Zl, Z2], etc., and 
references therein). However, up to now, to the best of our knowledge, the study of phase- 
transitional viscoelasticity has been carried out only for phenomenological models [SI, S2, 
S3, SS, Z2] for one-dimensional shear flow, with classical double well potential and artificial 
viscosity-capillarity terms. The treatment of the general, physical, case was cited in [BLeZ], 
Appendix B, as an important direction for further study. 

In this paper, continuing the work of Antman and Malek-Madani [AM], Slemrod [SI, 
S2, S3], Schecter and Shearer [SS], and Barker, Lewicka and Zumbrun [BLeZ], we study the 
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existence of planar viscoelastic periodic traveling waves, compressible or incompressible, for 
deformations of arbitrary dimension, starting from the most general form of the physical 
equations. It turns out for general viscoelasticity model with surface energy terms that, 
similarly as observed for the phenomenological models studied previously, the periodic trav- 
eling waves can only be standing waves and the corresponding ordinary differential equation 
(ODE) system exhibits Hamiltonian structure. 

As we know, for a planar Hamiltonian system, we can use phase-plane analysis to study 
its closed orbits. This corresponds to the case of one-dimensional deformations, for which 
the unknowns are scalar in the ODE system (4.5) that we obtain. For higher dimensional 
Hamiltonian systems, this method does not apply directly. In order to prove the existence 
of non-constant periodic waves when the unknowns are vectors, we consider the problem 
under the framework of calculus of variation. There are several problems to overcome. 
First, we need to formulate the problem in proper Banach spaces. It turns out that the 
proper space for this is the periodic Sobolev space with mean zero property. Working in 
this framework amounts to prescribing the mean of the unknown over one periodic (no real 
restriction, since each periodic wave has a mean as long as it exists). 

Second, we need to make sure that the waves we find are not constant waves. We 
overcome this issue by considering the equations satisfied by the difference between the 
original unknown and its mean. This makes the element in our working space always a 
critical point, which helps to eliminate the possibility that the solutions we find are trivial. 

Third, in the global model of viscoelasticity, we need to consider the physical assumption 
(see Section 2 or 3), namely T3 > 0. This makes our variational problem a constrained one 
(see Remark 6.3). 

This kind of constraint usually leads to a variational inequality and is related to an 
obstacle problem. Meanwhile, this constraint makes our admissible set (to make the wave 
physically meaningful) not weakly closed. However, the asymptotic behavior of the elastic 
potential function (see (2.4)) will help overcome the related problem. We note that in the 
pure elastic case without surface energy and viscosity, this restriction on T3 imposes signifi- 
cant challenge in the mathematical analysis (see the discussions in [Ba, AN]). However, the 
additional regularity imposed by the surface energy term removes this dificulty. 

Besides the Hamitonian structure of the standing wave equations, we prove that for 
the general physical model, there exist non-constant periodic waves no matter whether the 
unknowns are scalar or not (see Theorem 6.12) under assumptions on the mean vector and 
periodic of the wave. For some specific phase-transitional models, we given explicit condi- 
tions under which the non-constant oscillatory waves exist (see Section 7). In particular, 
for the one dimensional models, we use phase-plane analysis to get detailed information on 
the wave phenomena (existence of periodic, homoclinic, heteroclinic waves, see Section 8). 
The issue of their stability falls into the framework of [HZ, JZ] because the systems are 
effectively strictly parabolic (in the sense that it may be transformed to a strictly parabolic 
system; see Appendix B) and hence are completely understood. 

Comparing our results with others (see the references of this paper), the problems here 
is interesting enough even only from the modeling point of view, without even finding any 
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waves. In [AM], the authors treated the shear flow without the surface effect and with 
isotropic assumption preventing phase transition model (see discussion in [BLeZ]). Here we 
consider models with surface energy and the materials are anisotropic, which gives rich wave 
phenomena. Antmann, Slemrod and others (see [A, AM, SI, S2, S3] and references therein) 
have previously studied phenomenological ID phase-transitional models with double-well 
potentials. Here we justify those types of qualitative models by direct derivation from the 
physical shear flow model; see Section 5.7. 

It is well known that there are many outstanding open problems in elasticity and vis- 
coelasticity ranging from steady state to dynamical ones (see [Ba]). Due to surface energy 
effect and associated higher regularity, we can give a neat treatment of our problem. It 
would be very interesting to explore by numerics cases that do not fit the hypotheses here 
(which are sufficient but by no means necessary in the non-scalar case) but nonetheless 
support periodic waves and also to explore either numerically or analytically the spectral 
stability of these waves. We hope to address these issues in a followup work [BYZ]. 

2 Viscoelasticity models with surface energy term 

In this section, we will follow the presentations of [AN, Ba, BLeZ, NPT]. Let be the 
reference configuration which models a viscous elastic body with constant temperature and 
density. A typical point in will be denoted by X. We use £ : n x R + ->• M 3 to denote 
the deformation(i.e., the deformed position of the material point X). Consequently, the 
deformation gradient is given by F := Vx£, which we regard as an element in M 3x3 . 

Adopting the notations above, the equations of isothermal viscoelasticity with surface 
energy effect are given through the following balance of linear momentum modified with 
surface energy 

(2.1) & - V x • (pW{V£) + -Z(V£, V6) - £ (V 2 0) = 0. 

We make the following physical constraint on the deformation gradient (see [Ba, BLeZ] 
and [AN, NPT] for the physical background), prohibiting local self- impingement of the 
material: 

(2.2) detF>0. 

In (2.1), the operator Vj- stands for the divergence of an approximate field. As in [D, NPT], 
for a matrix-valued vector field, we use the convention that the divergence is taken row-wise. 
In what follows, we shall also use the matrix norm \F\ = (tr(F T F)) 1 / 2 , which is induced 
by the inner product: F\ : F2 := tr(F^ F2). 

In view of the second law of thermodynamics (see [Ba, PB]), the Piola-Kirchhoff stress 
tensor DW : M 3x3 — > R 3x3 is expressed as the derivative of an elastic energy density 
W : E 3x3 ->• !+. Throughout the paper, we assume as in [AN, Ba, NPT] the elastic energy 
density function W which is the stored-energy function for the viscoelastic material is frame- 
indifference. Let SO(3) be the group of proper rotations in R 3 . Then the frame- indifference 
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assumption can be formulated as 

(2.3) W{RF) = W(F), VF£R 3x3 , Vi? € SO(3). 

Also, the material consistency (to avoid interpenetration of matter, (2.2), [AN, Ba]) 
requires the following important assumption: 

(2.4) W(F) -> +oo as det F -»• 0. 

We emphasize that viscous stress tensor Z : R 3x3 x R 3x3 — > R 3x3 depends on both the 
deformation gradient F and the velocity gradient Q = Ft = V£t = Vt> , where r = £$. From 
physical point of view, the stress tensor should also be compatible with principles of 
continuum mechanics: balance of angular momentum, frame invariance, and the Claussius- 
Duhem inequality. For the related mathematical description and corresponding stress forms 
see [AN, Ba, BLeZ]. 

The surface energy £ is given by 

for some convex density VP : R 3x3x3 — > R, compatible with frame indifference. 
The corresponding inviscid part of system (2.1) 

(2.5) £« - V x ■ (DW(V0) = 
can be written as 

3 

(2-6) (F,T) t + J2dx i (Gi(F,T))=0. 

i=i 

Above, (F, t) : Q — > R 12 represents conserved quantities, while Gi : R 12 — > R 12 are given by 

-Gi(F, t) = T x &i T 2 ei r 3 ei 

are the fluxes, and denotes the i-ih coordinate vector in R 3 . 
The convex density ^ contributes to equation (2.1) the term 

(2.7) V x • (f(V 2 0) = V x • {V x • DV(V 2 0}- 

In view of the orders of differentiation and convexity of ^, we may assume that 
$ > 0; *(0) = 0; D*(0) = 0; Sid < D 2 V(-) < Mid 



o i ° 

W(F) 



8F, 



1...3 



fc=i 
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where 5, M are two positive real numbers. The mapping relations (ignoring physical con- 
straints) are 



si : M 3x3x3 



^ . R 3x3x3 R 3x3x3 



DH : R 3x3x3 £(R 3x3x3 ; R 3x3x3 ) 

and Id is an element in the space £(M 3x3x3 ; R 3x3x3 ). When the operator Vx- reduces to 
the operatore d x where x is a one dimension variable, (2.7) takes the form 

(2.8) d x {d x D*{%Z)}. 
If we identify £ x as r, then d 2 £ = t x and (2.7) becomes: 

(2.9) d x {d x D*(r x )} = d x {D 2 V(T x )T xx }. 

Note that D 2 ^ : R 3 ->• C(R 3 ; R 3 ) when V x - reduces to X . So we assume that D 2 *(-) as 
matrix function satisfy the following assumption: 

(2.10) Sid < D 2 V(-) < Mid as operators. 

These mapping relations enable us to decide the proper form of surface energy term con- 
tributed to equation (2.1). Keeping this in mind, we study specific equations in the following 
sections. 



3 Equations and specific models 

In this paper, we focus on the interesting subclass of planar solutions, which are solutions in 
the full 3D space that depend only on a single coordinate direction; that is, we investigate 
deformations £ given by 

= X + U(z), X = (x,y,z), U = (U 1 ,U 2 ,U 3 )eR 3 . 

Corresponding to the above deformation or displacement £, the deformation gradient is 
given by taking the gradient with respect to X 

/i o U hz \ /l rA 
(3.1) F = 1 U 2 , z = 1 r 2 . 

\o 1 + U 3 J \0 t 3 J 

We shall denote V = (r, u) = (n, r 2 , r 3 , u\, u 2 , u 3 ), where t\ = Ui tZ , t 2 = C/ 2 , 2 , r 3 = 1 + C/ 3j2 
and ni = Ui t t,U2 = U 2 j,u 3 = ^3,t with the physical constraint 

T 3 >0, 

corresponding to det F > in the region of physical feasibility of V. 
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(I 1 ° Tl \\ 

Writing W(t) = Wl I 1 T2 J ) , we see that for all F as in (2.1) there holds 
\0 r 3 / y 

V x • (DW(F)) = (D t W(t)) z . 

That is, the planar equations inherit a vector-valued variational structure echoing the ma- 
trix valued variational structure (note that the left hand side is the divergence of DW(F)). 

In this paper, we first study the problems (traveling wave ODEs, Hamiltonian ODEs, 
existence of standing waves) for general elastic potential energy and give a rather general 
abstract existence result. Then we study local models by specifying the related terms in 
system (2.1) as follows: 

1. Viscous elastic potential W. A typical model, as described in [FP], has the form: 

W(F) = \F T F — C- | 2 • \F T F - C+| 2 , 

which is a potential for anisotropic material with material frame indifference property. We 
study models involving this phase-transitional elasticity potential(See Appendix B for re- 
lated computations). It is important to notice that the elastic potential here does not satisfy 
the asymptotic behavior when det F — > + . Especially this is an irrelevant assumption for 
shear models. Hence the related models are local models for the real physics. 

2. Viscous stress tensor Z. We use the following viscous tress tensor which is 
compatible with the principles of continuum mechanics (see [BLeZ]) 

Z(F,Q) = 2(det F)sym(QF- 1 )F- 1 ' T . 

We note that the related Cauchy stress tensor T2 = 2(det F)~ 1 ZF T = 2sym(QF~ 1 ) is the 
Lagrangian version of the stress tensor 2symVv written in the Eulerian coordinates. For 
incompressible fluids 2div(symVv), giving the usual parabolic viscous regularization of the 
fluid dynamics evolutionary system. 

3. The surface energy £. A typical example of the energy density function is 
*(P) = \\P\ 2 , so that 

£ (V 2 £) = V x • V 2 e = A X F, 

which is an extension of the ID case of [SI]. For simplicity, we take this surface energy 
to study the wave phenomena. We see that it is the surface energy term that makes the 
models have abundant wave phenomena. 

The system. In order to be consistent with the usual presentation, we use x G I 1 as 
the space variable instead of z, i.e. z — > x. With the above choices, we have the system 



(3.2) 



r t -u x = 0; 

u t + (t{t) x = {b(T)u x ) x - (d(T x )T xx ) x . 
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Here a = -D t W{t), d(-) = D 2 V(-), 

(I 0\ 

(3.3) d(-) = Id and 6(r) = T3" 1 1 . 

\0 2/ 

We are interested in the existence of periodic traveling waves of the above system. Note 
that our system (3.2) in the variable (r, u) 6 W 1 © R n (n = 1, 2, 3.) enjoys some interesting 
structures: it is a hyperbolic-parabolic system with partial parabolicity; it involves third 
order terms because of the surface energy effect. See Appendix B for discussions. 



4 Traveling-wave ODEs 

We seek a traveling wave solution of the system (3.2): (r(x,t),u(x,t)) := (T(x—si),u(x—st)) 
where s £ R is the wave speed. Let us denote in the following ' as differentiation with respect 
to x — st. For convenience, we still use x to represent x — st (Indeed, we will show that s = 
is necessary for the existence of periodic or homoclinic waves). With further investigation 
in mind, we write the related equations for the general class of viscoelastic models. Now 
from system (3.2), we have the ODE system 



(4.1) 



-st 1 - v! = 0; 

-sv! + <j(t)' = {b{r)u'Y - {d(T')T")'. 



Plugging the first equation into the second in the above system, we have the following 
second-order ODE in r: 

(4.2) sV + <t(t)' = -{b{r)sT')' - (d{r')T"y. 
In view of d(-) = D 2 ^(-), we readily see: 

(4.3) sV + u{t)' = -{b{r)sr'y - {DH{t')t")' . 
Choosing a specific space point, say xo, we integrate once to get: 

(4.4) s 2 t + <j{t) +q = -s6(t)t' - Df(r')' 

Here q is an integral constant vector. Relating this with the elastic potential function W, 
we have 

(4.5) -DW{t) + s 2 T + q = -s6(r)r' - Df(r')' 
Carefully note that the integral constant vector is given by 

(4.6) q = {DW{t) - s 2 t - s6(r)r - D*(t / )'} 



X=Xq 
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5 Hamiltonian systems 



Defining G(P) := (P,DV(P)) - *(P), we see that = (P,D 2 *). Here P £ R" and 
* : R n — >■ R is the surface energy functional (for our purpose in this paper, n = 1,2,3.), 
G : R n — > R a different scalar potential type function. Now we are ready to state a structural 
property about the traveling wave ODE system (4.5). 

Proposition 5.1. When s = 0, the system (4.5) is a Hamiltonian system with factor 
^P 2 ^^')^ , preserving the Hamiltonian integral 

H(t, t') = -W(t) + qr + G(r') = constant. 

Here W can be taken in particular as the phase-transitional elastic potential (see Appendix 
A and also p. 36 of [BLeZj), with (01,02,03) = (ti,T2,ts)) 

W(r) = ( 2 (r 2 + r| + e 2 ) + (|r| 2 - 1 - e 2 ) 2 )* - 16e 2 r 2 , 

and ty(p) as any convex function with ^(0) = and d^(0) = 0, simplest case ^(P) = 
|P| 2 /2. (Note: in this simple case G(P) = ^>(P)-) 

Proof. When s = 0, the traveling wave ODE (4.5) becomes: 
(5.1) -dW(r) + q= -DV(t')' 

and the constant q = {DW(t) — D^(t')'} . Noticing the positive-definiteness of D 2 ^(), 

X=XQ 

we may write the ODE as a first order system by regarding r, r' as independent variables: 

t' = t'= [D^irT'D^ir'y; 
[ ' ' t" = -[D 2 V(t j )]-\-DW(t) + q) 

Now, consider the energy surface given by: 
(5.3) H(t,t') := - W{T) + qT + G{T'). 

We see that 



(5.4) 



^-H(T,T') = -DW(r)+q. 
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Comparing (5.2), (5.4), we see that the traveling wave ODE is a Hamiltonian system 
with factor 7 := [D 2 ^(r / )] _1 . Thus, (4.5) preserves the Hamiltonian H. We can see this 
also by the explicit computation, writing Q = x — st: 

±H(T,T>) =^H(T,T>)T> + ^H(T,T>)T" 



dr v ' ' dr 1 



dr' 



8t' 



dr 



□ 



From the above structural information, we easily get a necessary condition for the exis- 
tence of periodic or homoclinic waves, extending results of [OZ] in a one-dimensional model 
case. 

Theorem 5.2. For (4.5) with s ^ 0, there holds dH/d( «s 0, where 

s 2 

(5.5) H(r, t') := -W(t) + — |r| 2 + qr + G(t'), 

so that no homoclinic or periodic orbits can occur unless s = 0. 

Proof. Considering the evolution of -^H(t, t') along the flow of traveling wave ODE system 
(4.5), we have 

= {-D t W(t) + q + s 2 r, t') + {DG(t'), t") 
= (-D t W(t) + q + s 2 r, t') + {D 2 ^(t')t',t") 
= (-D t W(t) + q + s 2 r, t') + (D 2 V(t')t", t') 
= (-D t W(t) + q + s 2 T + DV(t')', t') 
= (-s6(t)t',t'). 

The conclusion thus follows from the positive definiteness of 6(r). □ 

The Hamiltonian system. From the above analysis, we may set s = 0. The traveling 
wave ODE system reduces to the following form with an integral constant q 



(5.6) 



-t" = -D t W{t) + q; 
q = {D t W(t) - t"} 



X=X() 



If we take the Hamiltonian point of view, the corresponding Hamiltonian for the above 
system is 

H(t,t>) = ±\t>(x)\ 2 + V(t,t>), 
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where V(t,t') := q ■ t(x) — W(t(x)). The periodic solutions of the system are confined to 
the surface H(t,t') = constant. 

In the following, we list the elastic potential and related information for the phase- 
transitional models we deal with in this paper for completeness and future study. To get 
these models, we fix one or two directions of r as zero, or, in the incompressible case, T3 = 1 
(as described in [AM, BLeZ], the latter is an imposed constraint, that is compensated for 
in the T3 equation by a Lagrange multiplier corresponding to pressure. 

We refer the reader to [BLeZ], Section 3, for details of the derivations of these models. 



5.1 2D Incompressible Shear Model 

This model corresponds to setting T3 = 1. 

(5.7) W(t) = (2rl + 2(r 2 - e) 2 + (|r| 2 - e 2 ) 2 ) (2r 2 + 2(r 2 + ef + (|r| 2 - e 2 ) 2 ) . 
Its gradient components are 

(5.8) D T1 W(r) = 8n(|r| 2 + 1 - £ 2 ){2(|r| 2 + e 2 ) + (|r| 2 - e 2 ) 2 }; 

(5.9) D T2 W(t) = 8t 2 (\t\ 2 + 1 - e 2 ){2(|r| 2 + e 2 ) + (|r| 2 - e 2 ) 2 } - 32r 2 e 2 . 

The Hessian components are 

(5.10) 

w u := D T1T1 W(r) = 8(|r| 2 + 1 - e 2 + 2r 2 ){2(|r| 2 + e 2 ) + (|r| 2 - e 2 ) 2 } + 32r 2 (|r| 2 + 1 - e 2 ) 2 ; 

(5.11) w 12 = W21 ■= D t1T2 W(t) = 16T 1 T 2 {2(|r| 2 +e 2 ) + (|r| 2 -e 2 ) 2 } + 32r 1 T 2 (|r| 2 + l-e 2 ) 2 ; 
(5.12) 

w 2 2 := D T2T2 W(t) = 8(|T| 2 + l-e 2 +2r 2 )[2(|r| 2 +e 2 ) + (|T| 2 -e 2 ) 2 ] + 32[T 2 (|r| 2 + l-e 2 ) 2 -e 2 ]. 

5.2 ID Shear Model I: r 3 = 1; r 2 = 0. 

The elastic potential becomes 

(5.13) W(t)= (2r 2 + 2e 2 + (T 2 -e 2 ) 2 ) 2 . 
The first order derivative is 

(5.14) D T1 W(r) = 8n(r 2 + 1 - e 2 ){2(r 2 + e 2 ) + (r 2 - e 2 ) 2 }. 
The second order derivative is 

(5.15) D T1T1 W(r) = 8(3rf + 1 - e 2 ){2( Tl 2 + e 2 ) + (rf - e 2 ) 2 } + 32r 1 2 (r 1 2 + 1 - e 2 ) 2 . 
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5.3 ID Shear Model II: r 3 = 1; r x = 0. 

Correspondingly, the elastic potential becomes 

(5.16) W(r) = (2(75, - ef + (r| - e 2 ) 2 ) x ( 2 (r 2 + e) 2 + (r| - e 2 ) 2 ) 
The first order derivative is 

(5.17) D T2 W{t) = 8r 2 (r 2 + 1 - e 2 ){2(r| + e 2 ) + (r| - e 2 ) 2 } - 32r 2 e 2 ; 
The second order derivative is 

(5.18) D T2T2 W(t) = 8(3r 2 + 1 - e 2 ){2(r 2 + a 2 ) + (r| - e 2 ) 2 } + 32{r 2 (r 2 + 1 - e 2 ) 2 - e 2 }. 

Remark 5.3 (ID vs. 2D shear solutions). Evidently, solutions of 1-D shear models I or 
II determine solutions of the full 2 shear model obtained by adjoining r 2 = or t\ = 
respectively. It is worth noting that the structure dW{r) = r/(|r|) + c(0, r 2 ) T for f a scalar- 
valued function and c a scalar constant yields that the only solutions t(x) of (5.6) with 
i] ■ t = c 2 = constant for some constant vector n are those satisfying c 2 /(|r|) + c?7 2 r 2 = 
constant, which gives by direct computation r = constant or else rj • r = and ?] 2 t 2 = 0, in 
which case t\ = or r 2 = 0. That is, the ID systems derived here are the only solutions 
of the 2D shear model that are not genuinely two-dimensional in the sense that they are 
confined to a line in the T-plane. In particular, if the mean of T\ or r 2 over one period is 
not zero, then we can be sure that the solution is genuinely two-dimensional. 

5.4 ID Compressible Model III 

In this case t\ = r 2 = and we denote r = T3. The potential and its derivatives are given 
below. The elastic potential becomes 

(5.19) W(t) = (2£ 2 + (t 2 -1-£ 2 ) 2 ) 2 
The first order derivative is 

(5.20) D T3 W(t) = 8t 3 (t 2 - l - e 2 ){2e 2 + (r| - l - e 2 ) 2 }; 
And the second order derivative is 

(5.21) W33 := D T3T ,W(t) = 8(3rf - 1 - e 2 ){2e 2 + (r| - 1 - e 2 ) 2 } + 32t 2 (t 2 - 1 - e 2 ) 2 . 



5 HAMILTONIAN SYSTEMS 



12 



5.5 Another two 2D compressible models. 

First, we consider the case r = (t 2 ,T3) t 6 R+. The elastic potential W and derivatives are 
as follows. 

(5.22) W(t) = (2(75 - e) 2 + (|r| 2 - 1 - e 2 ) 2 ) (2(r 2 + e) 2 + (|r| 2 - 1 - e 2 ) 2 ) ; 
The gradient components are 

(5.23) D T2 W(t) = 8t 2 (\t\ 2 - e 2 ){2(r| + e 2 ) + (|r| 2 - 1 - e 2 ) 2 } - 32r 2 e 2 . 

(5.24) D T3 W(t) = 8t 3 (\t\ 2 - 1 - e 2 ){2(r 2 + e 2 ) + (|r| 2 - 1 - e 2 ) 2 }; 

Similarly, we have the Hessian components 
(5.25) 

w 22 := D T2T2 W(t) = 8(\t\ 2 -e 2 + 2rf )[2(rf + e 2 ) + (|r| 2 - 1 - e 2 ) 2 ] + 32[r 2 (|r| 2 - e 2 ) 2 - e 2 ]. 
(5.26) 

^23 = ^32 := D T3T2 W(t) = 16T 2 T3{2(T 2 2 +£ 2 ) + (|T| 2 -l-e 2 ) 2 }+32T 2 T 3 (|r| 2 -e 2 )(|r| 2 -l-e 2 ) 
(5.27) 

^33 := J D T 3 T 3^(r) = 8(|T| 2 -l-e 2 + 2r 2 ){2(r 2 + e 2 ) + (|r| 2 -l-e 2 ) 2 } + 32r 2 (|T| 2 -l-e 2 ) 2 

Second, if we fix the r 2 direction and let r := (t±, T3) , we get another 2D compressible 
model. We omit the details here as the form is obvious. 

5.6 The full 3D model. 

In this case r = (ri, r 2 , t^) t G K+. corresponding to the phase-transitional elastic potential 
function W, we list the components of D 2 W(t) := (u>r,)3 X 3- 

(5.28) w u = 8(\r\ 2 + 2r 2 - £ 2 ){2(|r| 2 - r 2 + e 2 ) + (|r| 2 - 1 - e 2 ) 2 } + 32r 2 (|r| 2 - e 2 ) 2 

(5.29) wia = W21 = 16rir 2 {2(|r| 2 - r 3 2 + e 2 ) + (|r| 2 - 1 - e 2 ) 2 } + 32tit 2 (|t| 2 - e 2 ) 2 

(5.30) W13 = w 31 = 16r 1 r3{2(|T| 2 -T3 2 +e 2 ) + (|T| 2 -l-e 2 ) 2 }+32T 1 r 3 (|T| 2 -e 2 )(|r| 2 -l-£ 2 ) 

(5.31) w 23 = w 32 = 16r 2 r3{2(|T| 2 -T3 2 +e 2 ) + (|T| 2 -l-e 2 ) 2 }+32T 2 r 3 (|T| 2 -e 2 )(|r| 2 -l-£ 2 ) 

(5.32) w 33 = 8(\r\ 2 + 2r| - 1 - e 2 ){2(|r| 2 - r| + e 2 ) + (|r| 2 - 1 - e 2 ) 2 } + 32r 2 (|r| 2 - 1 - 6 2 ) 2 

(5.33) „ 22 = 8(|r| 2 - e 2 + 2r 2 )[2(r 2 + r 2 + e 2 ) + (|r| 2 - 1 - e 2 ) 2 ] + 32[r 2 (|r| 2 - e 2 ) 2 - e 2 }. 
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Remark 5.4. Similarly as in Remark 5.3, we find that solutions of the 2D compressible 
models are genuinely two-dimensional, in the sense that they are not confined to a line in 
the r-plane, unless they are solutions of the ID model derived above: in particular, the mean 
over one period of (ti, T2) is zero. Likewise, solutions of the full 3D compressible model are 
genuinely three-dimensional in the sense that they are not confined to a plane, unless they 
are solutions of one of the 2D models derived above, in particular, the mean of t\ or of T2 
over one period is zero. 

5.7 Justification of phenomenological models 

We note that, for the case of 1-D shear flow, the viscosity coefficient given by (3.3) become 
b,d = constant, and the elastic potential W is of a generalized double- well form. Thus, we re- 
cover from first principles the type of phenomenological model studied in [SI, S2, S3, SS, Z2], 
though with a slightly modified potential refining the quartic double- well potential assumed 
in the phenomenological models. The 2-D shear flow gives a natural extension to multi- 
dimensional deformations, which is also interesting from the pure Calculus of Variations 
point of view (see the following section), as a physically relevant example of a vectorial 
"real Ginzberg-Landau" problem of the type studied on abstract grounds by many au- 
thors. Finally, we note that the various compressible models give a different extension of 
the phenomenological models, to the case of "real" or nonconstant viscosity. 

6 Calculus of Variations. 

In this section, we formulate the problem in the framework of Calculus of Variations and 
give the proof of the existence result. 

6.1 Space structure 

As a first step, we recall the notions of Sobolev spaces involving periodicity and introduce 
the space structure we are going to use (see [MW]). For fixed real number T > 0, let Cj? 
be the space of infinitely differentiable T-periodic functions from R to W 1 (for our purpose 



n= 1,2,3.). 



Lemma 6.1. Let u,v £ L 



(0,T;R n ). If the following holds: for every f € CI 




then 
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and there exists a constant vector c in 1 N such that 

u(t) = / v(s)ds + c a.e. on[0,T]. 
Jo 

The function v := u' is called the weak derivative of u. Consequently, we have 

u(t) = [ u'(l)dl + c, 
Jo 

which implies the following: 

u(0) = u(T) = c; 

u(t) = u(s) + u'{l)dl. 

Proof. For the mean zero property, we could consider the specific test function f = ej. For 
the integral formulation, we could consider the use of Fubini Theorem and Fourier expansion 
of / to conclude ([MW]). □ 

Define the Hilbert space H\ as usual (hence reflexive Banach space) with this inner 
product and corresponding norm: for 

(u, v) := / (u,v) + (u',v')ds; 
Jo 

\u\\ 2 := / \u\ 2 + \u'\ 2 ds. 
Jo 



Next, we collect some facts for later use. 
Proposition 6.1. (Compact Sobolev Imbedding property) TVf CC C[0,T] compactly. 
Proposition 6.2. If u G H], and (1/T) J Q T u(t) dt = 0, then we have Wirtinger's inequality 

[ T \u(t)\ 2 dt < (T 2 /4tt 2 ) [ T \u'(t)\ 2 dt 
Jo Jo 

and a Sobolev inequality 

\u\l<(T/12) [ T \u'(t)\ 2 dt. 
Jo 

The Compact Sobolev imbedding property will give us the required weak lower semi- 
continuity property for the nonlinear functionals. The Wirtinger's inequality supplies us 
equivalent norms in related Sobolev spaces with mean zero property (see [MW] for complete 
proofs). 
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6.2 Variational formulation of the problems 

Now given real positive constant T, we consider problem (5.6) in Hj, 

f-r" = -D t W(t) + q = -D t {W{t) -q-r); 
\r(0) - r(T) = 0;r'(0) - r'(T) = 0. 

Let us first consider the cases and formulations without the physical restriction T3 > 0. 
Assume that: 



1 r 

t := — r{x)dx = m. 
T Jo 



Here m £ M n , n = 1,2,3 and we will use bar to represent mean over one period similarly. 
Hence, we consider the following problem 



(6.1) 



V'(x) = DW(t) - q; 
< t(0) = t(T);t>(0) = t>(T); 
,tIo r(x)dx = m. 



If we seek periodic solutions, q can be determined by integrating the equations above 
over one period; that is, 



1 t 1 

q = - I DW{r(x))dx. 
T Jo 



1 T 

Define v(x) = t{x) — m. We see easily that y J v(x) dx = 0, and v(x) satisfies the 
system of equations: 

v"(x) = DW(v + m) — q; 

(6.2) < v(0) = v(T);v'(0) = v'(T); 

,tSqv(x) dx = 0. 

For later convenience, we rewrite the above system as 

'v"{x) = DW(v + m) - DW(m) + DW(m) - q; 

(6.3) lv(0) = v(Ty,v'(0) = v'(Ty, 

^f T v(x)dx = 0. 

Define W(v) = W(v + m) — DW{m) ■ v and q = q — DW(m). We get the following 
problem 



( „ji 



v"{x) = DW(v) - q; 
(6.4) I v(Q) = v(T);v'(0) =v'(T); 

^J T v(x)dx = 0. 

Here q can be determined by integration: ^ Jq DW(v) dx = q. 
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Remark 6.3. Imposing the physical restriction t% > 0, we require that the T3- component 
of the mean vector m be positive: 777.3 > 0. Hence we need the assumption v% > —777,3 f or 
models involving T3 direction on [0,T]. 

Define F(v) = W(v + m) — W(m) — DW(m) ■ v and introduce the functional 

rT -1 pT 



T(v)= [ -\v'\ 2 dx+ [ F(v)dx 
Jo 2 Jo 



'0 

on the space 

H^ :={v€H^,v=^J vdx = 0}. 

Proposition 6.4. is always a critical point of the functional! defined above on H^ . 
Proof. It is easy to verify that for <f> G H^ Q , there holds 

I'(v)<p= [ v' ■ <f>' + DW(v) ■ <pdx. 
Jo 

Taking v = and noticing that DW(0) = 0, we get the desired result. □ 

Remark 6.5. By our formulation, we make always a critical point and it corresponds to 
the constant solution. This geometric property supplies us nice way to exclude the possibility 
that the periodic solution we find is constant, i.e., to help prove the periodic waves we find 
are oscillatory. 

Proposition 6.6. Without physical restriction on T3, the critical point of I corresponds to 
the solution of (5.6). 

Proof. This can be regarded as a simple consequence of Corollary 1.1 in [MW], For com- 
pleteness, we write the details here. First, assume v solves 

V'(x) = DW(v) - q; 
(6.5) I v{0) = v{T);v'{0) = v'{T); 

,Tl T v(x)dx = 0. 

Multiplying the equation by 4> G H\ and integrating to get 

f v'<t> +DW(v)-(j)dx = 0, 
Jo 

i.e v is a critical point of X. 

Next, we assume v is a critical point and <f> 6 IVf. Then (j) — cj) £ H-ro- Hence we have 

f v' ■ (</) - 4>)' + f DW{v) ■ {(f> - 4>) dx = 
Jo Jo 
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i.e. 



f v'-4>'+[ DW(v)-</>- [ DW{v) ■ (j)dx = 
Jo Jo Jo 

Noting that 4> = ^ f (j)dx, we find that the left-hand side expression above is: 

[ v'-cf)'+[ DW(v) ■(/)-[ DW(v)-(^ f c/)dx)dx = 
Jo Jo Jo T J 

Noticing that ^ J Q DW(v) dx = q, we get 

rp rp 

[ v' ■</>'+[ (DW(v)-q)-(f>dx = 0, 
Jo Jo 

which implies v" = DW{v) — q. □ 

Remark 6.7. If we consider models involving the restriction V3 > —m^, we need to consider 
a variational problem with this constraint. 

In view of the physical properties of the elastic potential function W (in particular, the 
polynomial structure of the phase-transitional potential function), we can apply the direct 
method of the calculus of variation to show existence. In order to deal with the integral 
constant q, we may restrict the admissible sets (or choose proper function space) on which 
we consider the functional or use Lagrange multiplier to recover it by adding restriction 
functional on the original space on which the functional is defined. 



In the following, we state with proofs some propositions of general nonlinear functionals. 
These propositions and further materials can be found in [De, Ni, ZFC] and the references 
therein. We will use them to prove the existence of periodic waves. 

Proposition 6.8. Let X be a Banach space, I a real functional defined on X and U be a 
sequentially weakly compact set in X. If I is weakly lower semi- continuous, then I attains 
its minimum on U, i.e. there is xq £ U, such that I(xo) = inf^gjy I(x). 

Proof. Let c := inf xg [/ I{x). By definition of inf, there exists {x n } C U such that I{x n ) — > c. 
In view that U is sequentially weakly compact, {x n } admits a weakly convergent subse- 
quence, still denoted by {x n }. Denote xo G X the corresponding weak limit. Since U is 
weakly closed, we know xo £ U. Noticing that weakly lower semi-continuity of /, we have 
c = lim n /(x n ) > I(xq). By the definition of c, we in turn know I(xo) = c > —00, which 
completes the proof. □ 

It is well-known that a bounded weakly closed set in a reflexive Banach space is weakly 
compact. In particular, a bounded closed convex set in reflexive Banach space is weakly 
compact since weakly close and close in norm are equivalent for convex sets. Hence we have 
the following corollaries: 
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Corollary 6.9. Let U be a bounded weakly closed set in a reflexive Banach space X and 
I be a weakly lower semi- continuous real functional on X . Then there exists xq G U such 
that I(xo) = inf^gj/I. 

Definition 6.10. A real functional I on a Banach space X is said to be coercive if 

lim I(x) = +oo. 

|x|*-»-+oo 

Corollary 6.11. Any coercive weakly lower semi- continuous real functional I defined on a 
reflexive Banach space X admits a global minimizer. 

6.3 A general existence result 

In this part, we first give a general result for models with the physical assumption r 3 > 0, 
i.e, v 3 > — m 3 . We will assume the following conditions on the potential W 

(Al) W G C 2 and W{t) -»• +00 as r 3 -»• 0+. For r 3 < 0, define W{t) = +00; 
(A2) There exist a positive real number C such that W(t) > —C for r G R n (n = 
1,2,3); 

(A3) There exists a constant vector m G := {m G K 3 ; m 3 > 0} such that a{D 2 W(m)}n 
Ml / 0. Here cr{D 2 iy(m)} is the spectrum set of D 2 W(m). 

Because of our formulation, we only need that the potential is bounded from below. For 
a potential such an assumption is quite natural. Assumption (A3) amounts to saying that 
there is a point where the potential is concave. From the physical point of view, this is 
quite reasonable. 

Theorem 6.12. // (^) 2 < X(m), then we have a physical nonconstant periodic wave 
solution for the problem (3.2) for which the mean over one period of r is m. Here —X(m) 
is the smallest eigenvalue of D 2 W(m). 

Define two subsets of LL\ by 

A\ := {v £ H^ ;v 3 > -m 3 }; 

A 2 ■= {v G Ht j0 ;v 3 > -m 3 }. 
Remark 6.13. The admissible set A\ is not weakly closed in H^ . 
Lemma 6.14. Under assumptions (A1)-(A3), T is a coercive functional on H^ Q . 

Proof. By the definition of 1, we just need to consider the part Jq F(v) dx. By assumption 
(A2), we have 
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/ F(v)dx= \ W(v + m) — W(m) — DW(m) ■ v dx 
Jo Jo 

= / W(v + m) — W(m) dx 
Jo 

> [ -C-W(m)dx 
Jo 



= - (C + W(m))T > -oo. 

□ 

By the above lemma, we see for sufficient large R the minimizers of X on Ai are restricted 
to the sets Ai := Ai n B H i^ [0, R] for i = 1, 2 where -B^a [0, JS] is the closed ball with center 

and radius R in H^ . Define Si := {v G Ai\T(v) = inf^g^ Obviously, we have 
5i := {f G A;I(w) = inf eeA X(5)}. 

Lemma 6.15. ^ is a weakly compact set in H^ . 

Proof. A2 is bounded by its definition. Since H\ is reflexive, we know A2 is weakly se- 
quentially compact. Also, A2 is convex. Indeed, we can use the definition of convexity of a 
set to check this easily. An appeal to Sobolev embedding theorem yields that A2 is closed 
in norm topology of Hf. For a convex set, closeness in norm topology and weak topology 
coincides, hence we have that A2 is weakly closed. Putting this information together, we 
have shown that A2 is weakly compact. □ 

Lemma 6.16. Z is a weakly lower semi- continuous functional on H^ Q . 

Proof. Let v n — > v weakly in H^ . By Sobolev imbedding, we have v n — > v uniformly 
in [0, T}. Hence we have J Q T F(v n )dx — > F(v) dx. Because of the mean zero property, 
Iq \v'\ 2 dx is of norm form, hence it is a weakly lower semi-continuous functional. □ 

Lemma 6.17. 5*2/0 and V3 > —m^ + e for v G S2 under the assumption of Theorem 6.12. 
Here e is a positive constant. 

Proof. By Proposition 6.8, S 2 / 0. Note that G A2, 1(0) = and hence Z(v) < 0. Hence 
we will have V3 > —m^ + e. Or else by Sobolev imbedding there would be a positive mea- 
surable set of [0,T] on which T3 > is vanishingly small. By assumption on the asymptotic 
behavior of W with respect to T3 — > 0, we would have I(v) = f^ (l/2)\v'\ 2 dx + J Q T W(v + 
m) — W(m) dx > = X(0), a contradiction. □ 



Lemma 6.18. ^ Si = S2 under the assumption of Theorem 6.12. 
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Proof. Consider the second variation. An easy computation shows that for v,(j> in H^ Q 

l"(v):((t)®(f>)=[ \<j)'\ 2 dx+ I' D 2 W{v + m):{(f)®<j))dx. 
Jo Jo 

To show ^ S2, consider 

X"(0) : (0® <j>) = [ \(f)'\ 2 dx+ [ D 2 W(m) : (<f> ® </>) dx . 
Jo Jo 

Let <p(x) = r/sin(^^) for < 77 < 7723 and vq £ I 3 be a unit eigenvector corresponding to 
— A(m). We see that 4>{x) := 4>(x)vq £ Ai- Since is a critical point of X on H^ and 

l"(0):(f ®W= [ T v 2 & 2 (cos(^)) 2 dx- A(m) f T r] 2 (sm(^)) 2 dx 



Jo T T jo ± 

=^{(y) 2 " A (-)}<°- 

Hence we see that S2 and Si = S2 is obvious. □ 

Proof of Theorem 5.11. Combining Lemma 6.14-Lemma 6.18, we finish the proof of Therem 
6.12. □ 

Remark 6.19. The condition (^") 2 < A(m) in Theorem 6.12 on the period T, is readily 
seen by Fourier analysis to be the sharp criterion for stability of the constant solution 
t = m, u = 0. Equivalently, it is the Hopf bifurcation condition as period is increased, 
marking the minimum period of bifurcating periodic waves. Thus, it is natural, and no 
real restriction. On the other hand, there may well exist minimizers at whose mean m 
W is convex; this condition is sufficient but certainly not necessary. Likewise, there exist 
saddle-point solutions not detected by the direct approach. 



6.4 Relation to standard results, and directions for further study 

In the scalar case tGR 1 , the condition that D 2 W{m) have a negative eigenvalue is equiva- 
lent to convexity of the Hamiltonian H at the equilibrium (m, 0), under which assumption 
there are many results on existence of periodic solutions of all amplitudes; see, for example, 
[Ra] and later elaborations. Likewise in the vectorial case r 6 M. d , d > 1, if D 2 W(m) < 0, 
then we may appeal to standard theory to obtain existence of periodic solutions by a variety 
of means; indeed, the convexity condition may be substantially relaxed for solutions in the 
large, as described in [Ra], and replaced by global conditions ensuring, roughly, star-shaped 
level sets of the Hamiltonian. On the other hand, review of the potentials considered here 
reveals that, typically, it is a single eigenvalue of D 2 W that becomes negative and not all 
eigenvalues, and so these methods cannot be directly applied. 

It is an interesting question to what extent such standard methods could be adapted 
to the situation of a Hamiltonian potential (in our case — W) with a single convex mode. 
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Existence of small amplitude periodic waves at least is treatable by Hamiltonian Hopf 
bifurcation analysis. The question is to what extent if any one can make global conclusions 
beyond what we have done here, in particular, to relax for large solutions the nonconvexity 
condition on W at m. Finally, it would be interesting to find natural and readily verifiable 
conditions for existence of saddle-point solutions in this context. 

7 Existence of periodic solutions for specific models 

In this section, we focus on the existence of periodic waves for the incompressible models 
(i.e. models with T3 = 1), namely the ID models I, II and the 2D incompressible model. 
First note that these models do not involve the T3 direction. Hence we have no condition 
corresponding to (Al). However, the specific phase-transitional elastic potential energy 
function W has good growth rate when |r| — >■ +00 for r = t\,T2 or (ti,T2). This will make 
our functionals coercive. Hence we have the following: 

Theorem 7.1. For the incompressible models, there exist non-constant periodic standing 
waves respectively if the mean m (either vector or scalar) satisfies {AS) and (^) 2 < X(m). 
When m is scalar, assumption (A3) means m lies in the elliptic region of the viscoelasticity 
system (2.1). 

Proof. It is easy to see the corresponding functionals are coercive, weakly lower semi- 
continuous functionals on the reflexive Banach spaces H^ - Hence Corollary 6.11 applies. 
The verification that the global minimizers respectively are not zero is entirely the same as 
in Lemma 6.18 by considering the second variation. □ 

Next, we specify the corresponding conditions in Theorem 7.1 for these incompressible 
models. 

7.1 ID Shear Model I 

The condition is 

(7.1) (^) 2 < -8(3m 2 + 1 - e 2 ){2(m 2 + e 2 ) + (m 2 - e 2 ) 2 } - 32m 2 (m 2 + 1 - e 2 ) 2 . 
In particular, if m = 0, the condition reads 

(7.2) (^) 2 <-8(l-e 2 )(2 £ 2 + £ 4 ). 

Condition (7.2) illustrate that our assumption is not a void assumption. Also, in the mean 
zero case, (7.2) holds only if e > 1. Comparing this with the existence result by phase-plane 
analysis (see in particular section 8.1), we see that these results match very well. 
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7.2 ID Shear Model II 

The condition is 

(7.3) (^) 2 < -8(3m 2 + 1 - e 2 ){2(m 2 + e 2 ) + (m 2 - e 2 ) 2 } - 32{m 2 (m 2 + 1 - e 2 ) 2 - e 2 }. 
In particular, if m = 0, the condition reads 

(7.4) (^) 2 < -8(1 - e 2 )(2e 2 + e 4 ) + 32e 2 = 8(e 6 + e 4 + 2e 2 ). 

Condition (7.4) implies in particular that for any e > 0, we have long-periodic oscillatory 
waves. Similarly, for any given T > 0, we have oscillatory waves as long as e > large 
enough. Comparing with the phase-plane analysis (see section 8.2), we see that the related 
obtained wave phenomena match very well. 

7.3 2D Incompressible Shear Model 

In this case D 2 W(m) is given by its components 

Wu ■= D nn W(m) = 8(|m| 2 + l-e 2 +2m 2 ){2(|m| 2 +e 2 ) + (|m| 2 -e 2 ) 2 }+32m 2 (|m| 2 + l-e 2 ) 2 ; 

Wl2 = W21 : = D TlT2 W(m) = I6m 1 m 2 {2(\m\ 2 + e 2 ) + (\m\ 2 -e 2 ) 2 } + 32m 1 m 2 (\m\ 2 + l-e 2 ) 2 ; 

w 2 2 ■= D T2T2 W{m) = 8(|m| 2 +l-e 2 +2m 2 )[2(|m| 2 +e 2 )+(|m| 2 -e 2 ) 2 ]+32[m 2 (|m| 2 +l-e 2 ) 2 -e 2 ]; 

The corresponding condition is (^") 2 < A(m). This is obviously a rather mild condition. 
To see this, we can consider in particular the mean m = (mi,0) T or (0, m2) ■ Then the 
results on the two ID incompressible models readily give the conclusion because we have 
diagonal matrices. 

Based on the analysis of these conditions, we have in particular (m = case): 

Theorem 7.2. For the 2D shear model, ID shear model I and II, we have the following 
existence result of periodic viscous traveling/standing waves: 

(1) Given any e > eq > 0, for any T satisfying T > T(e) > 0, system (5.6) hence (2.1) 
has a nonconstant periodic solution with some appropriate integral constant q; For the ID 
model II we have e$ = 0. 

(2) Given any T > 0, for any e satisfying e > e(T) > 0, system (5.6) hence (2.1) has a 
nonconstant periodic solution with some appropriate integral constant q. 

Remark 7.3. From the above theorem, we see in particular indeed for the 2D shear model, 
we have infinitely many nontrivial periodic viscous traveling waves with appropriate corre- 
sponding q values. In particular, we have a sequence of waves with minimum positive period 
T -»• +oo. 
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Remark 7.4. Remarks 5.3 and 5.4 show that solutions of the specific 3D model of Section 
5 with m±, mi, and 771,3 nonzero are genuinely 3- dimensional in the sense that they are not 
confined to a plane in r-space, and that solutions of the various 2D models of Section 5 
are genuinely two-dimensional if the means of both components are nonzero. That is, we 
have constructed by the variational approach solutions that are not obtainable by the planar 
phase-portrait analysis of the ID case (see just below). On the other hand, a dimensional 
count reveals that, generically, the periodic solutions nearby a ID solution are all ID, and 
likewise the periodic solutions nearby a 2D solution are all 2D. 

8 ID existence by phase-plane analysis. 

In this section, we discuss how to generate periodic waves for ID models. In (5.6), the 
integral constant q = D t W(t-) — r'L- Here (t_,t") is the vector evaluated at some specific 
space value xo- If there indeed exist periodic-T waves, q = ^ J Q T DW{t{x)) dx. 

By the variational formulation and the usual bootstrap argument, we conclude that 
the periodic waves are classical solutions of the system (5.6). For the possible homo- 
clinic/heteroclinic solutions of (5.6), standard elliptic regularity theory also guarantees that 
these waves are classical solutions. There must be points (say xq,x\) in a period [0,T] 
such that t'(xo) = and t"(x\) = , etc., if such periodic solution did exist for r scalar. 
The reason is that r(x) cannot be always monotone and convex in view of periodicity (this 
applies to all derivatives). Hence we can make the integration constant have the form 
q = D t W{t(x\)) for convenience a priori. Then we can show existence, which in turn guar- 
antees the a priori assumption. Hence, we could assume q = D t W(t-) to show existence. 
We adopt this convention in the following analysis. 

The guiding principle is that the ODE systems are planar Hamiltonian systems. To get 
complete and clear pictures of the phase-portraits, we just need to specify the "potential 
energy" term V(t,T-) in the Hamiltonian H(t,t'). 

8.1 ID Shear Model I: r 2 = 0; r 3 = 1. 

In this section, we denote r = t\. We use similar notation in other sections. Recall the 
elastic potential W(t) = I 2r 2 + 2e 2 + (t 2 — e 2 ) 2 J and its first and second order derivatives 



D T1 W(r) = 8n(r 2 + 1 - e 2 ){2(r 2 + e 2 ) + (r 2 - e 2 ) 2 }; 
D t1ti W(t) = 8(3rf + 1 - e 2 ){2(r 2 + e 2 ) + (r 2 - e 2 ) 2 } + 32r 2 (r 1 2 + 1 - e 2 ) 2 . 



The traveling wave ODE and corresponding Hamiltonian system are 



(8.1) 



r " = W'(t) - W'(t-). 



(8.2) 
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We write the Hamitonian system as follows: 



(8.3) 



T 



H(t,t') - V(t;t- 



) = E-V{t,t-). 



2 



Here E are constants corresponding to energy level curves of H(t,t') and V(t,t-) := 
qr- W(t). 

First, we determine the number of equilibria of the Hamiltonian system, hence focus on 
the solution of W'(r) = q. 

Note that W'(r) is an odd function on the real line, hence we just need to study its 
graph on the interval (0, oo). In view of the expression of W (r), we need consider the cases: 



For the case < e < 1, we have W"(r) > for r real, hence W'(r) is strictly monotone 
increasing and 



Hence in this case, for any given r_, the solution of W'(t) = W'(r_) is r_ and unique. 
Similar analysis holds true for e = 1. Considering the definition of V(t,T-), we have 

Proposition 8.1. When < e < 1, V(t;t-) has exactly one critical point r_, which must 
be a global maximum. 

Remark 8.2. In this case, our Hamiltonian system admits no periodic orbit for any r_ (or 
equivalently, for any q). 

Next, consider the case e > 1. In this case, we can see from the expression of W'(t) 
that W'(t) has three distinct zeros: —\Je 2 — 1, 0, \J e 1 — 1. A qualitative graph of W (t\) 
is as follows: 



(1) < e < 1; (2) e > 1. 



W'(0) = 0; W'(r) >0, forr>0; W'(r) < 0, forr<0; 
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The corresponding graph for the potential V(j\,T\ -) is as follows 




Proposition 8.3. The function W'(t) has exactly two critical points. 

Proof. By symmetry, we do the following computations: Denote r\ := X and e 2 := a > 1. 
We want to show that the function 



f{X) := (3X + 1- a)[2(X + a) + (X - a) 2 ] + AX(X + 1 - a) 2 

has exactly one zero when X > 0. 

First, noting that /(0) < and /(pp) > 0, we know that f(X) has a root on (0, pp). 
Also note that /(X) > on [pp,oo), hence we just need to show that f{X) admits a 
unique zero on (0, pp). Computing the derivative, we have 



f'{X) = 3[7X 2 + 10(1 - a)X + a 2 + 2a + 2(1 - a) 2 ]. 

Denote 

A = 100(1 - a) 2 - 28[a 2 + 2a + 2(1 - a) 2 ). 

If A < 0, we know that f'{X) > 0, hence f(X) is monotone increasing, which implies 
that f(X) admits a unique zero; 

If A > 0, we will have two positive roots for f'(X) = and the smaller one is 1Q ( a ~^~^ _ 

However, we can show that 10( - a ~ 1 1 ^~ v/ ^' > ^p, hence the function f(X) is monotone in- 
creasing on the interval (0, pp), which also implies the uniqueness of the zero. 

Now we have a clear picture on the potential W'(r) (see the graph for W'{t\)). 



□ 



Proposition 8.4. When e > 1, the function W'(t) is an odd function with 3 zeros and 2 
critical points and goes to infinity when r — > +oo. 
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Denote the two critical values of W'(t) as q* > and — q*, for convenience denoting 
Q = q*. Then we have the following property: 

Proposition 8.5. Assume e > 1. When \q\ > Q, the equation W'(t) = q has exactly one 
solution; When \q\ = Q, the equation W'(t) = q has exactly 2 solutions; When \q\ < Q, the 
equation W'{r) = q has exactly 3 solutions. 

As the solutions of W'(t) = q correspond to the critical points of V(t;t-), we have: 

Theorem 8.6. For \q\ > Q, the Hamiltonian system admits no periodic orbit; For \q\ < Q, 
the Hamiltonian system admits a family of nontrivial periodic orbits. Further if q = 0, the 
Hamiltonian system also admits a heteroclinic orbit. 

Proof. For \q\ > Q, we know that V(t;t_) has a global maximum and hence the Hamilto- 
nian admits no periodic orbit. For the case \q\ < Q, we know that the potential V(t;t_) 
must have exactly 3 critical points with 2 local maxima and 1 local minimum. Also, we 
know that F(r;r_) has strictly lower energy at the local minimum than at the two local 
maxima. Hence the existence of a family of periodic orbits follows. Further, when q = 0, 
the energies at the two local maxima of V(t; t_) are the same, hence we get an heteroclinic 
orbit. □ 

Remark 8.7. We may compare the two energy values o/y(r;r_) at the two local maxima. 
If they are equal (when q = in particular), we have a heteroclinic orbit. In general, they 
are not equal to each other, which yields a homoclinic orbit. 

8.2 ID Shear Model II: n = 0; r 3 = 1. 

Recall the elastic potential 



D T2 W(t) = 8r 2 (r| + 1 - e 2 ){2(r| + e 2 ) + (r| - e 2 ) 2 } - 32r 2 e 2 ; 
D T2T2 W(t) = 8(3rf + 1 - e 2 ){2(r 2 + e 2 ) + (r| - e 2 ) 2 } + 32{rf (rf + 1 - e 2 ) 2 - e 2 }. 



As above, we list the Hamiltonian system and the potential V(t;t-). The system is: 



W(t) = (2(r 2 



e) 2 + (T 2 2 -e 2 ) 2 )x(2(T 2 + e) 2 + (T 2 -e 2 ) 2 ) 



and the relevant derivatives 



(8.4) 




The potential is V(r; r_) = qr — W(r). 
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Remark 8.8. There is a slight difference with the ID shear model I in the function W'(t). 
Because of this difference, we do not need to restrict the positive number e to get periodic 
orbits for the parameter q in proper range. The conclusions are completely the same when 
1 > e > as in ID shear model I when e > 1. 

We have the following: 

Proposition 8.9. When 1 > e > 0, the behavior of the function W'(r) is the same as that 
of the function W'(t) in the ID shear model I when e > 1. In fact, W"(t) is monotone 
increasing in this case for r > 0. 

Remark 8.10. For the range e > 1, numerics suggest that the behaviors are also the same 
as we may show that the function W'(t) has exactly three solutions and two critical points. 
We have a small problem to verify this by direct computation though we just need to show 
that f{X) > evaluated at the larger root of f'(X) (f(X) is defined similar as in ID shear 
model I as the second derivatives of the two potentials differ with a constant 32e 2 . Even 
without this, we still can conclude the existence of periodic orbits since the potential V(t; t_) 
admits a minimum. Together with the existence obtained by variational argument, we know 
that there are still infinitely many nontrivial periodic waves for any e > 0. 

8.3 ID Compressible Model III 

For this model, we need to pay special attention to the physical restriction T3 > when we 
do the phase-plane analysis. To find physical waves, we use a continuity argument and a 
simple comparison criterion. 

In this case t\ = T2 = 0; let r = T3. The potential and its derivatives are given below. 

The elastic potential becomes W(t) = ^2e 2 + (r| — 1 — £ 2 ) 2 ^ . Its first and second order 

derivative are D ta W{t) = 8T 3 (r 3 2 - 1 - e 2 ){2e 2 + (r| - 1 - e 2 ) 2 } and 

^33 := D T3T:i W(r) = 8(3rf - 1 - e 2 ){2e 2 + (r| - 1 - e 2 ) 2 } + 32t 2 (t 2 - 1 - e 2 ) 2 . 

We write V = V(r,q,e) := qr — W(t) in this section to emphasize the analytical 
dependence of V on the parameters q and e (because V is a polynomial). As in previous 
sections, we see: 

Proposition 8.11. (1)W'{t) = q always has one, two or three roots when \q\ > Q, \q\ = Q 

or \q\ < Q for some positive Q. In the case that W'(r) = q has 3 distinct roots, we denote 
them from small to large by ti , r m and r r . 
(2) W'(t) has exactly two critical points. 

As before, in order to analyze the existence of periodic or homoclinic/heteroclinic waves, 
we just need to consider the potential energy V(t, q, e). Further, in order to have physical 
waves, we need necessarily that — Q < q < 0. In this situation, the two roots r m and r r of 
W'(t) = q are positive. Noticing that r = r m is a local minimizer of V(r,q,e), there is a 
periodic annulus around r m . Hence we have the following proposition: 
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Proposition 8.12. When —Q < q < 0, there always exists a periodic annulus. 

To show existence of physical homoclinic orbit, we just need to compare the values of 
V(0,q,e) =: V(0) and V(j r ,q,e) =: V(r). We have 

Proposition 8.13. Let —Q < q < 0. IfV(0) > V(r), there is a physical homoclinic orbit; 
IfV(0) < V(r), there is no physical homoclinic orbit. 

In particular, for the case q = 0, the 3 distinct roots of W'{t) = q are easily seen to be 
n = -\/TT7Vm = and r r = -n. So y(0,0,e) = -(2e 2 + (1 + e 2 ) 2 ) 2 < V(r r ,0,e) = 
— (2s 2 ) 2 . By continuity, we have the following conclusion: 

Proposition 8.14. There exists a constant r\ > such that if—r] < q < 0, then there exists 
no physical homoclinic orbit. 

Proof. When q = 0, V(0,0,e) < V(r r ,0,£). Thus, by continuous dependence and Proposi- 
tions 8.11 and 8.13, we have the relation V(0) < V(m) holds when q < is small and the 
conclusion holds. □ 

Next, we study the existence of physical homoclinic waves when — Q < q < is large. 
For this purpose, we first set e = then proceed by a perturbation argument. When e = 0, 
the corresponding elastic energy function and its derivatives are: 

W(t) = (r 2 - l) 4 ; 

W'(t) = 8r(r 2 - l) 3 ; 
Ty"(r) = 8(T 2 -l) 2 (7r 2 -l). 
Note that W"(t) = has roots r = ±1, ±y/l/7 (this can be easily computed). 

For this potential V(r,q,e = 0) := qr - W(t), we need \q\ < W'(-^/l/7) = 8(^ f^ 

to have a homoclinic wave. For physical ones, we need — 8(j) 3 ^J^ < q < 0. Consider the 

case q — > — 8(|) 3 -^/^ from the right, we see the largest root T m of W'(t) = q tends to \f\Jl . 
Consequently, V(t = 0,q, e = 0) — > — W(0) = — 1 and the right local maximum value of 
V{T m ,q,e = 0) -»• -8(f f^\-W{^Jlj7) < -1. By Proposition 8.11 and Proposition 8.13, 
we have: 

Proposition 8.15. For the compressible ID model, assume —Q < q < (Q as in Propo- 
sition 8.11). Then, when e > and q + Q are small, we have a physical homoclinic orbit. 



9 Time- evolutionary stability 

We conclude by discussing briefly the question of time-evolutionary stability of viscoelastic 
traveling waves, up to now not considered. 
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9.1 Spectral vs. nonlinear stability 

A very useful observation regarding the earlier phenomenological models r t — u x = 0, u t + 
dW(r) x = bu xx — dr xxx , b, d > constant, for 1-D shear flow, made by Schecter and 
Shearer in [SS], was that for a wide range of b, d, specifically, d < 6 2 /4, the system can be 
transformed by the change of independent variable u — >■ u := u — ct x , c{b — c) = d to the 
fully parabolic system 

T~t U x — CT XX , 

u t + dW(r) x = (b - c)u xx , 

thus allowing the treatment of nonlinear stability by standard parabolic techniques, taking 
into account, for example, sectorial structure, parabolic smoothing, etc. 

Quite recently, this observation has been profoundly generalized by M. Kotschote [K], 
who showed that a somewhat different transformation in similar spirit may be used to 
convert viscoelastic or fluid-dynamical equations with surface energy (resp. capillarity) 
effects to quasilinear fully parabolic form, in complete generality, not only to the cases 
4d > b 2 previously uncovered for the phenomenological model, but to the entire class of 
physical models considered here. For further discussion/description of this transformation, 
see Appendix B. 

This reduces the question of nonlinear stability to a standard format already well studied. 
In particular, it follows that (except possibly in nongeneric boundary cases of neutrally 
stable spectrum) nonlinear stability is equivalent to spectral stability, appropriately defined. 
This follows for heteroclinic and homoclinic waves by the analysis of [HZ], and for periodic 
waves by the analysis of [JZ]. 1 For precise definitions of the notions of spectral stability, we 
refer the reader to those references; in the shock wave (heteroclinic or homoclinic) case, see 
also the discussion of [BLeZ]. Spectral stability may be efficiently determined numerically 
by Evans function techniques, as in for example [BHRZ, BLeZ, BHZ, BJNRZ1, BJNRZ2]. 
We intend to carry out such a numerical study in a followup work [BYZ]. 

9.2 Variational vs. time-evolutionary stability 

More fundamentally, perhaps, there is relation between variational stability of periodic 
waves and their time-evolutionary stability as solutions of (2.1). In particular, the energy 
functional that we minimized in constructing periodic solutions is essentially the self-same 
functional that defines the mechanical energy of the system, a Lyapunov functional that 
decreases with the flow of (2.1); for further discussion, see Appendix C. This gives a strong 
link between the two notions of stability. Indeed, it can be used to directly show that the 
periodic waves constructed as minimizers of the associated variational problem are time- 
evolutionarily stable with respect to co-periodic perturbations (Remark C.2, Appendix C). 

Moreover, as discussed further in [BYZ], the indirect spectral arguments of [Z2] on 1- 
D heteroclinic and homoclinic waves extend to to the general-dimensional and or periodic 

lr The analysis of [JZ] concerns modulational stability, or stability with respect to localized perturbations 
on the whole line; co-periodic stability may be treated by standard semigroup techniques [He]. Spectral 
analyses of [OZ, BYZ, PSZ] suggest that modulational stability occurs rarely if ever for viscoelastic waves. 
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case, yielding the much stronger result that variational stability in each of these contexts is 
necessary and sufficient for time-evolutionary stability (co-periodic stability, in the case of 
periodic waves, and variational stability constrained by a prescribed mean). Moreover, these 
arguments yield at the same time the curious fact that unstable spectra of the linearized 
operator about the wave must, if it exists, be real. These properties give additional insight, 
and additional avenues by which time-evolutionary stability may be studied. 

In particular, this shows that the waves we have constructed are the (co-periodically) 
stable ones. However, these are not necessarily the only stable waves, as we did not con- 
struct all minimizers of the variational problem, but only those with mean satisfying a 
nonconvexity condition. Just recently (in particular, after the completion of the analysis of 
this paper), there has been introduced in [PSZ] a different, more direct argument showing 
equivalence of variational and time-evolutionary stability, which yields at the same time 
concise conditions for variational stability. These yield in particular that the sharp condi- 
tion for stability is not the condition of nonconvexity of W at m, defined as the mean over 
one period of r, but rather the "averaged" condition of nonconvexity of the Jacobian with 
respect to m of the mean over one period of DW(t). See [PSZ] for further details. 



A Appendix: Phase-transitional elasticity 



In this appendix, we collect some computations for the phase-transitional elasticity ([BLeZ],[FP] 
and references therein). 2 



W(F) :-- 



F F — C- 



F F — C 4 



where 



'I 
















1 












C± = {F T F)± := ( 1 ±e 

,0 ±e l + e 2 r 

/0 

Evidently, W is minimized among planar deformations at the two values A± = ±e 

\o o o 

Indeed, we have then 

W(F) = (2rl + 2(r 2 - e) 2 + (|r| 2 - 1 - e 2 ) 2 ) (2r 2 + 2(r 2 + e) 2 + (|r| 2 - 1 - e 2 ) 2 ) , 



'0 Tl 



1 Tl" 



where, as in [BLeZ], A = r 2 and F = 1 r 2 , so that 



,0 r 3 , 



,0 T 3 , 



(A.l) 



1 Tl 

F 1 F = | 1 r 2 

<n r 2 |r| 2 



2 For the convenience of the reader, we note that the the vector (ai, 02,03) in [BLeZ] corresponds to 
(ti,t 2 ,t 3 ) here. 
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and 
(A.2) 



W{F) = 




n 

\t\ 2 -1-s 2 



2 









I 









r 2 - e 



T 2 -e 

M 2 -i-e 2 



T 



(A.3) 



(F T F-C±f 



* (t"2 T e) 2 * 

* * rf + (r 2 Te) 2 + (|r| 2 



l-£ 2 ) 2 , 

B Appendix: The structure of the viscoelastic model. 

In this section, we verify the structural properties of the viscoelastic model and related 
modified systems obtained by the apporach of Kotschote [K]. To be clear, we collect these 
related systems. The original system is 



(B.l) 



Tt ~ U x = 

u t + (t(t) x = Q)(t)u x ) x - (d(T x )T xx ) x . 



Here a = -D t W{t), d(-) = D 2 $>(-), d(-) = Id and 
(B.2) 



1 N 
6(r) = r, 1 [ 1 

2; 



Introducing the phase variable z := t x , we may write as a quasilinear second-order 
system 



(B.3) 



Tt ~l" %x U x — T xx 

u t + a(r) x = (6(r)n a 



(d(z)z a 



Remark B.l. This transformation, introduced in [Kj, is similar in spirit to but more 
general than the one 3 introduced by Slemrod [SI, S2, S3] and used in [OZ] for an artificial 
viscosity /capillarity model. 

We can slightly modify the above system in the second equation. Then we have the 
following system 



(B.4) 



T~t Z x U x — T xx 
%t ~\~ %x — U xx + T xx 

u t + ct(t) x = (b(T)u x ) x - (d(z)z x ) x . 



A transformation (r, u) — > (r, u — ct x ) reducing the model to a parabolic system of the same size. 
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If we write the above system (C.3) and (C.4) in matrix form 

U t + f(U) x = (B(U)U x ) x 

using the variable U := \ z , the corresponding matrix B becomes: 




-d(z) b(r) 

Also, the corresponding matrix Df(U) for the two systems are: 




Proposition B.2. (Strict parabolicity) Systems (B.3) and (BA) are both strictly parabolic 
systems in the sense that the spectrum of B have positive real parts. 

Proof. Comparing the two matrices B above, we know they have the same spectrum. We 
prove this proposition for t,z,u € M 3 . The lower dimension cases becomes easier and the 
computations are totally the same . Pick one of the B's, say 

'h 3 3 

3 3 h 

K 3 -d(z) 6(r), 

To compute the spectrum of B, consider the characteristic polynomial det(Aig — B) = 0, 
which is 

/(A-l)/ 3 3 3 \ 

det 3 AI 3 -I 3 = 0. 

V 3 +d( Z ) X-b{T)J 

Doing Laplace expansion and elementary column transformation , we get 

A 2A 

det{(A - l)/ 3 } det{(A 2 / 3 - Afe(r) + d(z))} = (A - 1) 3 (A 2 - - + 1) 2 (A 2 - — + 1) = 0. 

T3 T3 
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Prom the first factor of the above degree 9 polynomial, we get three equal root 1 which 
has positive real parts . The other 6 roots also have positive real parts noticing that 
r 3 > 0. □ 

Proposition B.3. (Nonzero characteristic speeds) The corresponding first order systems 
of (C.4) has nonzero characteristic speed at r where the matrix D 2 W(t) are strictly positive 
definite. 

Proof. Again, we prove this for t,z,u G M 3 . To prove the corresponding first order system 
is non-characteristic, we consider the spectrum of the matrix I O3 J3 O3 I . We get 



the following system for the characteristic speed: 



K Da(r) 3 3 



A/ 3 — / 3 J 3 

det(A/ 9 - Df(U)) =| 3 (A - l)h 3 |=0. 

-Dg{t) 3 A/ 3 , 

By direct computation, we know: 

(B.5) det(A/ 9 - Df(U)) = det ((A - l)/ 3 ) det ( - Da(r) - A 2 / 3 ) = 0. 

It is easy to see we have three roots 1 which is not 0. The other roots satisfy the algebraic 
equation: det ^ - Do(t) - \ 2 I^j = det (^D 2 W(t) - A 2 / 3 ^ = 0. Hence the proposition 
follows. □ 

Proposition B.4. (Same spectrum) For system (BA) and its first order system, the matrix 
Df(U) and B~ l Df(U) have the same spectrum. 

Proof. We prove this for the variables r, z, u G M 3 . It is easy to verify that 

/ 3 3 3 
B- 1 = ( -b(r) b(r) -h 
-h h 3 



Since 



we immediately get 



3 h -h 
Df(U) = I 3 h 3 
k Do-(t) 3 3 

3 h ~h 
B^DfOJ) = I -Da{r) 3 6(r) 

3 3 h 
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Considering the corresponding eigenvalue problem, we have: 

. ( AI 3 -h h 

det (A/9 - B- l DfOJ) \ = det Da{r) XI 3 -b(r) 

\ 3 3 (A-l)/ 3/ 

Doing Laplace expansion and performing basic transformation, we get: 

det (A/3 - / 3 ) det (a 2 / 3 + £><t(t)) = det (\I 3 - J 3 ) det (a 2 / 3 - D 2 W{t)^ = 0, 
which implies the conclusion by noticing (C.5). □ 

C The principle of least action and relative entropy 

In this section, we discuss further the relation between variational stability and time- 
evolutionary stability with respect to co-periodic perturbations of periodic waves. It is 
well-known that one of the fundamental laws in nature is the principle of least action (also 
known as Hamilton's principle, Maupertuis' principle). The most general formulation of 
the law governing the motion of mechanical system can be derived from this principle. 

Consider a system which occupies at instants t\ and £2 positions given by the generalized 
co-ordinates vectors q x ,q 2 6 R N ■ Then the principle says that the system moves between 
such positions that makes the action functional S := J t * 2 L(q(t),q'(t),t) dt takes the least 
possible value. Here the function L in terms of the time t, the generalized co-ordinates q 
and velocities q'(t) is called the Lagrangian of the system. 

Based on the principle of least action, one can give a unified formulation of laws of mo- 
tion in mechanics in terms of differential equations (for example, wave equation, membrane 
equation, small oscillations, equations in fluid dynamics, etc). In particular, one can derive 
the Hamilton-Jacob theory of mechanics from which one could obtain systems of conser- 
vation laws (See [Lan]). Also, this principle is the foundation of calculus of variation and 
critical point theory (see [MW]). For more related details on this topic see [Lan, MW, PB] 
and the references therein. 

To be self-contained, we recall the system (3.2) (which is (C.l) below now) 



(C.l) 



T t -u x = 0; 

u t + ct(t) x = (b(T)u x ) x - (d(T x )r x 



Here the functions b, d are the same as before. From [BLeZ] we know the following is 
the associated mechanical entropy 

u 2 

(C.2) v(t,u) = — + W(t) + *(t x ) 



Consider the following mechanical energy for a given positive period T 
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r 1 r T u 2 

(C.3) £ = / 7](T,u)dx= / — + W(t) + *(r x ) dx 

In any case, the mechanical energy £ decreases along the flow for periodic boundary 
condition, with dissipation as follows 



f 

Jo 



d f 

—£{t,u) = / (uu t + DW{r)T t + DV(T x )T xt )dx 
at Jo 
rT 

-u x b(r)u x + u x d(r x )T xx + u x u(t) + DW(t)u x + D^(r x )u xx )dx 

io 

-T 

u x b(r)u x dx < 0. 

'o 

Hence we see immediately that periodic traveling wave solutions must have u = constant vector, 
from which we then find easily the speed s = in view of the relation u t = —su x = and 

T t ~U x = 0. 

Next, we adopt the periodic Sobolev space framework to discuss the relation between 
the least action functional we defined and the relative entropy 

(C.4) £(t, u;< t >, < u >) := £(t, u) - £(< r >, < u >) - D£(< r >, < u >) • (r, u), 



where < r >:= y J T r(x) dx as in physics literature and similarly for 
After a brief computation, we get that the relative entropy is given by 

£(t, u;< t >,<u>) = [ *(t x ) + W(t) - DW(< r >) • r - W(< r >) dx 

Jo 

f T !, ,2 1, ,2 

+ / ol u l — — | < ii > | — < u> -uax. 
io 2 2 

Prom the system (3.2) we know that the structure is preserved under the transformation 
u — > u + c where c is an arbitrary constant vector. Hence without loss of generality, we can 
let u = 0. Hence we get the following expression by further choosing < u >= m as before 

S(t, 0; 0, m) = [ V(t x ) + W(t) - DW(m) ■ r - W(m) dx 
Jo 

= [ -\t x \ 2 + W{t) - DW(m) -t -W{m)dx 
Jo 2 

Defining the translated variable v(x) = t(x) — m, we get the relation between the least 
action functional and the relative entropy 



u. 



£(v, 0; 0, m) = T(v) + constant, 
where the constant is given by T(DW(m) ■ m). 
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Remark C.l. The relative entropy is a rather common construction, meant to be stationary 
about the reference configuration (in our case < r >= m and < u >= 0). 

Remark C.2. The discussion above also sheds some light on the relation between time- 
evolutionary properties and variational structure of the Hamiltonian structure of our prob- 
lem (See [Z2J for further discussions). By the method of [GSSI, GSSII], we see that the 
mean- constraint minimizers we constructed are necessarily stable in the time- evolutionary 
sense with respect to co-periodic perturbation if they are stable in the variational sense. 
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